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Suppression of superconductivity due to non-perturbative saddle points in the 

nonlinear cr-model 

D. A. Pesin and A. V. Andreev 
Department of Physics, University of Washington, Seattle, WA 98195, USA 

We study superconductivity suppression due to thermal fluctuations in disordered wires using the 
replica nonlinear cr-model [NLgM). We show that in addition to the thermal phase slips there 
is another type of fluctuations that result in a finite resistivity. These fiuctuations are described 
by saddle points in NLaM and cannot be treated within the Ginzburg-Landau approach. The 
contribution of such fluctuations to the wire resistivity is evaluated with exponential accuracy. The 
magnetoresistance associated with this contribution is negative. 

PACS numbers: 74.78.Na, 73.21. Hb, 74.25.Fy, 74.40.+k, 74.25.Ha 
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Suppression of superconductivity in low dimensional 
samples has been studied for a long time. In particular, 
in quasi-one-dimensional samples supercurrent decay at 
temperature, T, below and near the critical temperature, 
Tc, is believed to occur via thermally activated phase 
slips (TAPS). This idea was put forward by Little [ij, 
and the quantitative theory was developed by Langer, 
Ambegaokar, McCumber and Halperin (LAMH) ^ yj. 
Within the LAMH theory the wire resistance below Tc, 
R, has activated behavior, R ex exp(— Ai^/T), where AF 
is the activation energy for a TAPS. 

Most experimental results appear to be in good agree- 
ment with the LAMH theory [J, |5j . However some dis- 
crepancies have also been reported |a,t3l- In particular, 
the negative magnetoresistance observed in thin Pb wires 
in Ref. |g] has no explanation within the LAMH theory. 

From a theoretical viewpoint, the LAMH theory relies 
on the Ginzburg-Landau description of thermal fluctua- 
tions resulting in finite resistivity. In this paper we study 
superconductivity suppression in disordered wires using a 
more general description of a superconductor in terms of 
the replica nonlinear cr-model (NLaM) [g] that includes 
both the pairing and the Coulomb interaction. This for- 
malism describes the low energy physics in terms of the 
order parameter. A, the fluctuating electric potential, 
V, and the Q-matrix that describes the fluctuating elec- 
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tron Green function. Since the Ginzburg-Landau func- 
tional can be obtained [9j from the NLaM the LAMH 
optimal fluctuation is automatically taken into account 
within this formalism and appears as a saddle point of 
the NLaM action. 

We show that due to the interplay of Coulomb inter- 
action and disorder there exist different saddle points 
of the NLaM action that provide alternative channels 
for supercurrent decay. These fluctuations cannot be de- 
scribed within the Ginzburg-Landau approach. They can 
be understood as follows. The interplay of disorder and 
Coulomb interaction leads to soliton-likc fluctuations of 
the Q-matrix analogous to those first studied in Ref. [iflj 
in the context of charge quantization in a single electron 
box with diffusive contacts. In a disordered wire these 
solitons are accompanied |ll| by a change of the imagi- 
nary part of the electric potential by 2tiT on the scale of 
thermal diffusion length, Lt — \/D/2'kT, where D is the 
diffusion constant. Below Tc the Coulomb and pairing in- 
teractions result in bound soliton-antisoliton pairs with 
a typical spatial separation, Lc ^ min{Ljv,L5}, with 
Ln - al^ln^^ and Lg - [p{Tc - T)AlAlT^]-\ 
where v is the density of states at the Fermi level, A is 
the wire cross section area, Aq is the equilibrium value 
of the superconducting gap, e is the electron charge, and 
d is the transverse dimension of the wire. In the vicinity 
of the critical temperature Lc ^ Lt, and the soliton and 
antisoliton within each pair are well separated. We show 
below that superconductivity is suppressed in the middle 
of the soliton-antisoliton pair, see Fig. ^ This results in 
an additional wire resistance 



5R{T) = Ro exp [-2G(Lt) ~ 2AF/T] 



(1) 



FIG. 1: Sketch of the electric potential, V(x), (solid line) and 
superconducting gap, A (2:), (dashed line) dependence on the 
coordinate x along the wire at the saddle point configuration. 



Here AF/T « 0.35G(LT)[(rc - T)/Tcf/^ is the acti- 
vation action of the LAMH phase slip, with G{Lt) = 
AttvDA/Lt ^ 1 being the dimcnsionlcss conductance of 
the wire segment of length Lt- Equation (^ is valid for 
G{Lt)-^/^ < (Tc - T)/Tc < G{Lt)-'^^^, where the first 
inequality corresponds to the Ginzburg criterion and the 
second ensures the condition Lc^ Lt- 

Although rigorous evaluation of the pre-exponential 



factor, _Ro: is rather difficult and left for future work, 
we speculate that the interior of each soliton-antisoliton 
pair acts as a segment of a normal metal, see discussion 
below Eq. UlUI) . This reasoning leads to the estimate 

Rq ~ Rn YJt , where Rjq is the wire resistance in the nor- 
mal state. 

The magnetoresistance associated with the contribu- 
tion of these saddle points is negative at low fields and is 
given by Eq. (|12|l . This may be relevant to the negative 
magnetoresistance observed in Ref. p . 

We consider a long disordered wire with many trans- 
verse channels. The temperature dependent coherence 
length ^(T) is assumed to exceed the transverse wire di- 
mension, d. The Thouless energy for the transverse mo- 
tion, Et = D/d?, satisfies the condition Ex ^ T. In 
this regime the wire is described by a one-dimensional 
NLaM. 

We write the replicated partition function for the sys- 
tem, {Z"') = (e^"~), with n being the number of repli- 
cas, as an imaginary time functional integral over the 
fermions. Upon disorder averaging, the Q- matrix (in 
the space of replicas and Matsubara frequencies) is in- 
troduced to decouple the disorder-induced quartic inter- 
action. Its matrix elements are 4x4 matrices in the 
space S (^ T, given by the product of spin, S, and time- 
reversal, T, spaces [lj,ll3. The auxihary fields Va{x,T) 
and Aa{x, r) decouple the Coulomb and pairing interac- 
tions in the replica a via Hubbard- Stratonovich transfor- 
mation. Proceeding in the usual way |3|, one arrives at 
the following form of the action for Q, T^ and A: 
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(2b) 



(2c) 



Sc = ]-ldTdxdx'Y,ya{x,T)C{x,x')Va{x',T).{2d) 
•' a 

Here Tr denotes the trace over the replica, Matsubara 
and S ^ T spaces. The matrices e, V and A have the 
following structure in the replica and S ® T spaces: e = 
tS'^^'Tsdr, V = J-^VoK, A = 5''\AaT+ - A^T,), whcrc 
T± — '^^^'"^^ and TiS are defined as t^ ^ ti® ua, with a^ 
and ti being the Pauli matrices in the S and T spaces re- 
spectively. The positive constant A denotes the attraction 
strength in the Cooper channel, and C{x,x') describes 
the inverse effective Coulomb interaction in the wire. In 
particular, for a homogeneous wire in the absence of a 
nearby gate its Fourier component is C{q) = 1/e^ In t5^. 
The saddle point equations are obtained by minimizing 
the action in Eq. (|2a|) with respect to A, y and Q. For 



this purpose one may neglect Sc in |5J if the number of 
channels in the wire is sufficiently large, e^vA 3> 1. This 
corresponds to the charge neutrality limit [ll| and gives: 



L»V(QVQ)-[£ + y + iA,Q] = 0, (3a) 

K- JtrQ^^(x)=0, (3b) 

A„+z^tr[r_g-(x)]=0, (3c) 

where tr is the trace in the S ®T space. Equation H3a|) 
is the Usadel equation, Eq. (j3b|) represents the charge 
neutrality condition, and Eq. l|Sc|) plays the role of the 
BCS self-consistency equation. 

Before discussing the superconducting case, let us 
briefly describe solutions of the saddle point equations Q 
in a normal metal, A^, = [ll|. Equations (|3a|l and (|3b|) 
possess a set of degenerate stationary spatially uniform 
solutions, where Q"^ = (5°''(5„™roSgn(e„ -|-27rrwa), K = 
2TTT'Wa with Wa being an integer, n, m being the Matsub- 
ara indices, and £„ — TrT{2n+l). The sum 4 ^^Wa = W 
(the factor 4 here arises from the 4x4 structure in 
the S (E) T space) defines the trace of the Q-matrix, 
TrQ = 2W. Different W^'s define topological classes of 
Q, i.e. Q matrices corresponding to the same W can be 
continuously rotated one into another in the replica and 
Matsubara spaces. Thus in addition to the above spa- 
tially uniform saddle point configurations, the action (j21) 
has spatially inhomogeneous finite-action solitons as sad- 
dle configurations. These solitons connect different Q's 
of the same topological class, and correspond to both 
Q and Va changing in space. These nonuniform saddle 
points are crucial for the subsequent considerations. 

Let us construct the soliton that connects Q"'' — 
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<5°^^nmToSgn(e„ — 2TiTwa) at X = oo, with w^'"^ = =f1, 
and all the other w's are zero. This corresponds to grad- 
ual change in the electric potential in replicas 1 and 2 
from zero at negative infinity to =p27rjr at positive infin- 
ity. For simplicity, we assume strong spin-orbit scattering 
and consider the Q-matrix belonging to the symplectic 
ensemble [l^l- In this case, the soliton Q-matrix can be 
parameterized by a single rotation angle 9{x) [lO|: 



Qn{x) 






(4) 



where u is the generator of rotation between the Mat- 
subara frequencies ttT in replica 1 and — ttT in replica 
2, uf^ = <5'^i<5^'<5„,o'5™,-i + <5'^2<5'i<5„,_i<5„,o. The in- 
dex i — Q,l labels diffuson-like, tq, and cooperon-like, 
Ti, rotations. Substitution of Eq. Q) into l|3b|) gives 
T4_2(2:) = T7i-T[l - cos6l(x)]. Then Eq. ^ gives 



v^e. 



1 



Zljrp 



sin 29^ = 0. 



(5) 



This equation has a solution 9i{x) = 2 arctan(e^/^^ 
which gives for the electric potentials 14.2(2;) = ^ttT[\ 



tanh(a;/I/T)]- The action (|2b|) of this sohton is G{Lt)- 
It can be shown that this saddle point is stable |l4|. 

In a macroscopic system, solutions with nonzero Wa's 
at X — > ±00 have infinite Coulomb action Ij2d|l and 
thus are forbidden. Inhomogeneous saddle points ap- 
pear as bound soliton-antisoliton pairs with Va{x) ap- 
preciably different from zero only in the interior of the 
pairs, see Fig. Q] In the regime of interest, G{Lt) 2> 1, 
the Coulomb action Ij2d|l does not affect the shape of 
the solitons, but merely determines the typical soliton- 
antisoliton separation, Ljv- The latter can be esti- 
mated by equating the Coulomb action of a pair, Sc '^ 
e-in{llh ' to unity, resulting in L^ ~ ^l^ In 5;^. We 
assume that Ln ^ Lt- This assumption corresponds 

to TTei < 21^(17)^ In^ -^^Wd ("^e' i'^ *h^ mean free time) 
and is always satisfied in the NLaM applicability region. 
Since the soliton and antisoliton in a typical pair are well 
separated, the pair action cost is 2G{Lt)- 

We now turn to the superconducting case. At T > Tc 
Eq. ij^ implies Aa — 0, and the solutions of Eqs. (|3a|l . 
(|3bl) remain the same as in the normal case. Below Tc, 
in the region of interest, Tc — T <C Tc, a small static 
superconducting order parameter A^ appears. The sad- 
dle point solution, Qs, can be expressed in terms of the 
normal state solution Q) as Qs — U^QnU, where U is 
a replica-diagonal rotation matrix in the Gorkov-Nambu 
space, U = S"-''U"-, with 



^ nm 



= <5„ . 



.cos^-T2e~'^«»<5„ 



-(m+i)Sgn(e„)sin^. 

(6) 
Here Xa are the order parameter phases. 

Solving Eq. \A&\ using the ansatz © with the spatially 
uniform Qjv — A, and substituting the result into Ij^cf) 
one recovers the Ginzburg-Landau equation for the order 
parameter, which contains both the homogeneous BCS 
solution and the LAMH saddle point. 

Let us turn to the inhomogeneous Qn, Eq. Q. In the 
replicas a 7^ 1, 2 the situation is the same as for Qm = A. 
For a = 1, 2 substitution of the ansatz © into Eqs. (O 
leads to the following equations for (j)'^, and A^, 



1 + cos 6* 2 
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sin (t)Q 2A 
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— cos 00, (7a) 
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n>0 



■ Sine 



(7b) 
(7c) 



Here we assumed that there is no current in the wire and 
set the order parameter phase, x, to zero. The additional 
equation for 6 aX Tc — T <^ Tc may be approximated by 
Eq. ((Sj. Therefore below we consider 6 in Eqs. (|7|) to be 
the (approximate) solution of Eq. Q corresponding to 
the soliton-antisoliton pair in the normal metal. With- 
out loss of generality we consider the soliton and anti- 
soliton to be positioned at a; = ±C with C ~ Lc ^ Lt, 



where Lc is the typical soliton-antisoliton separation in 
the presence of superconductivity, estimated below. 

The exact solution of the system of equations iQ is 
prohibitively difhcult. However, near the critical tem- 
perature the situation simplifies. Outside the soliton- 
antisoliton pair, \x\ > C, the GL approach may be used. 
From the viewpoint of the GL description the presence of 
the soliton-antisoliton pair imposes a boundary condition 
on the order parameter Aa at a; = ±C, 



VlnAa(2;)U=±C = ±k, 



(8) 



where n ^ L^ <C £,{T)^^. The form of this boundary 
condition can be obtained from the following considera- 
tions. The effective "critical temperature", T*, for the 
appearance of a static order parameter, Aq, is lowered 
inside the soliton-antisoliton pair. Indeed, for cos 6* = —1 
linearizing Eqs. I|7b|l and (|7c|l for a uniform static Aq we 
obtain 



27rz/A 



n>0 ' 



1 



(9) 



where as usual the sum over Matsubara frequencies in 
the r.h.s. must be cut off at the Debye frequency. The 
restriction n > in the summation arises because the 
last term in Eq. (f7c|l vanishes for cos 9 = —1. Expressing 
the l.h.s. of Eq. © in terms of Tc we find T^ = Tc/e^, 
where e ~ 2.71 is the base of the natural logarithm. Thus 
the regime of interest, Tc — T <^ Tc, corresponds to the 
"normal" phase inside the soliton-antisoliton pair [l5j. 
Since (T - T*)/T* - 1 the effective correlation length 
inside the pair is of order Lt, and the order parameter 
present outside the pair penetrates into its interior to dis- 
tances of order Lt ^ £,{T)- This leads to the boundary 
condition © similar to the proximity effect for a normal 
metal in contact with a superconductor jig . Then in the 
limit Lt ^ £,{T),h the boundary condition (jSJ may be 
replaced by Aa(±C) = 0, and the order parameter Aa{x) 
acquires the form. 



|Aa(x)|=Aoe(|x|-C)tanh 



C 



\/2C(T)' 



(10) 



where Q{x) is the step function. 

Since the order parameter profile outside the soliton- 
antisoliton pair coincides with that in the LAMH fluc- 
tuation in the zero current limit the action cost of such 
a configuration is Sq = 2G{Lt) -I- 2AF/T, where AF is 
the LAMH free energy barrier in the limit of zero bias 
current. The factor 2 in front of AF arises because the 
condensate suppression occurs in two replicas. In this 
estimate we neglected the Coulomb action and the con- 
densation free energy in the interior of the pair, 5'^ ~ 

eHnicfd) + CAv{Tc - T)Al/Tc. For a typical fluctuation 
S(^ ^ 1 <^ Sq. This gives for a typical fluctuation size, 
Lc ^ min{LAr,Ls}, where {Ls)~^ = y[Tc - T)AlA/Tl 



With exponential accuracy the contribution of such 
fluctuations to the wire resistance is determined by their 
action, which results in Eq. ^. Evaluation of the pre- 
exponential factor Rq is a difficult task. It may in prin- 
ciple be accomplished using linear response theory. This 
involves a cumbersome procedure of analytic continu- 
ation of the Matsubara susceptibility to real frequen- 
cies and is outside the scope of the present work. It 
is nevertheless plausible that the interior of the soliton- 
antisoliton pair acts as a segment of a normal metal |l7| , 
giving a contribution -^-(jh^jf- to the wire resistance. 
The number of soliton-antisoliton pairs in the wire can be 
estimated as -^^^e""^", where the prefactor corresponds 

to L/ Lt ways to accommodate a soliton, and Ld Lt pos- 
sibilities to put an antisoliton next to it. Following this 
argument one arrives at the estimate i?o — Rnjt-- 

We now discuss the magnetic field dependence of re- 
sistance Q . It arises from the decrease of condensation 
energy A_F due to magnetic field, and an increase of the 
action of cooperon-like solitons. In the presence of a mag- 
netic field, H, the NLaM action for a wire can be ob- 
tained [13 by changing {VQf -^ {V Qf - e^ {A^)[ts , QY 
in Eq. Ij2b|l . where (A^) is the square of the vector poten- 
tial, A, averaged over the wire cross-section (we assume 

that the magnetic length, Ih — a/ J77 j is much larger than 
the wire width, d). Each soliton or antisoliton constitut- 
ing a pair can be generated either by a diffuson-like or a 
cooperon-like rotation, z = 0, 1 in Eq. Q). The diffuson 
is insensitive to the magnetic field whereas the cooperon 
acquires a mass. As a result, Eq. ^ for the cooperon 
rotation angle on the soliton configuration becomes 



1 / 1 



Z \ IJrp 



1 



sin 20c = 0, 



(11) 



where L^ = e 



^'A^) ex H^ depends on the specific ge- 
ometry, in particular for a rectangular wire of width d, 
and a magnetic field perpendicular to it, L^ = Slj^/d^. 
The corresponding action change for a single cooperon- 
like soliton is G{Lt) (y/l + L'^/Ljj- iV The change 

in the condensation action, AF/T, may be evaluated us- 
ing the Ginzburg-Landau theory and is much smaller, of 
the order of G{Lt)^J{Tc - T)/T,. At relatively weak 

fields, L?fj > G{Lt)L'^ ^ L"^, it may be neglected while 
the action change for the cooperon-like soliton can be 
significant. In this regime we obtain |l8| 



6R{H) 



5R{T) 



\ _|_ g-G(LT)Lip/2L%, 



fl2) 



Here 5R[T) is given by Eq. (Q) and we took into account 
that there are four ways to assemble a soliton-antisoliton 
pair out of diffuson-like and cooperon-like solitons. From 



Eq. (|12|l we see that the resistance Q decreases with 
magnetic field. 

In summary, we have considered superconductivity 
suppression in thin wires due to non-pcrturbative sad- 
dle points in the NLaM . We showed that their presence 
leads to an additional (to the LAMH) contribution to 
the wire resistance below Tc and evaluated it with ex- 
ponential accuracy. This contribution is suppressed by 
a magnetic field, which may correspond to the negative 
magnetoresistance observed in Ref. J3|. 
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In principle, the soliton-antisoliton pair may be viewed 
as an SNS junction, but because its size, Lc exceeds the 
coherence length in the interior, which is of the order of 
Lt, the expected critical current is exponentially small. 
We neglect the magnetic field dependence of the pre- 
exponential factor. We believe this is justified for Lh 3> 

LT,aT). 



